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Abstract. Vanadium sesquioxide, V2O3, boasts a rich phase diagram whose 
description necessitates the accounting for many-body Coulomb correlations. Spectral 
properties of this compound have been successfully addressed within dynamical mean 
field theory to an extent that results of recent angle resolved photoemission experiments 
have been correctly predicted. While photoemission spectroscopy probes the occupied 
part of the one-particle spectrum, optical experiments measure transitions into empty 
states and thus provide complementary information. In this work, we focus on the 
optical properties of V2O3 in its paramagnetic phases by employing our recently 
developed "generalized Peierls approach" . Compared to experiments, we obtain results 
in overall satisfactory agreement. Further, we rationalize that the experimentally 
observed temperature dependence stems from the different coherence scales of the 
involved charge carriers. 



1. Introduction 

Vanadium sesquioxide, V2O3, has been the subject of extensive theoretical and 
experimental studies for now more than three decades. It is considered as the prototype 
compound, that undergoes a Mott-Hubbard transition [U [2] in its purest form. Indeed, 
the high-temperature (T>T]\g^g|) metal-insulator transition upon chemical substitution, 
(Vi_a;Cra;)203, is isostructural and no magnetic order is acquired. Early theoretical 
approaches resorted to the Hubbard model to explain the electronic properties of V2O3. 
However, over the years, experiments indicated that the physics of this material is more 
involved and a realistic multi-orbital setup is needed for the complexity of the correlation 
effects taking place. For reviews see e.g. [H [21 [3]. 

The field of correlated materials gained major momentum from the development 
of dynamical mean field theory (DMFT) Pj. In combination with standard density 
functional based methods like the local density approximation (LDA) the calculation 
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of spectral properties of materials with strong electronic Coulomb interactions became 
possible. Over the past years, LDA+DMFT enlightened our understanding of materials 
such as transition metals, their oxides or sulphides, as well as f-electron compounds [5]. 
Several works highlighted the applicability of the technique to V2O3 [6l [TJ [H [9l [TOl [TT] . 
In our previous work [H], we find that the metal-insulator transition is not due to the 
Brinkman-Rice mechanism [12] in its single-band form, but results from the impact 
of Coulomb correlations on the crystal-field splitting. Owing to its octahedral oxygen 
surrounding, the vanadium 3d-orbitals split into two and three lower lying t2g orbitals. 
The two manifolds of bands are isolated in energy both from each other and from other 
orbitals. The trigonal part of the crystal field further splits the t2g into an aig and 
two lower lying degenerate orbitals. The local Coloumb correlations result in an 
increased ai^-e^-splitting with respect to the LDA, causing a charge transfer that pushes 
aig spectral weight above the Fermi level. By computing momentum-resolved spectral 
functions [131 El E]) we made explicit predictions for angle-resolved photoemission 
experiments. Recent measurements on (Vi_a;Cri^)203 (x=0.011) [I5] nicely agree with 
the theoretical spectra, further validating our current understanding of this compound. 



2. Optical properties — prelude 

An experimental probe complementary to photoemission is optical spectroscopy, which 
is commonly analyzed in terms of the optical conductivity 

k 

that is given by a convolution of spectral functions 74k(a;)li|. The transitions are weighted 
by dipole matrix elements, called Fermi velocities, V]s_^a, and the Fermi functions f{ui) 
select the range of occupied and empty states, respectively. V is the unit-cell volume, 
a,/? denote cartesian coordinates, and 3fJ(T"^ is the response in a— direction for a light 
polarization along {3. Both spectral functions and velocities are matrices in the basis 
of localized orbitals, which we index by L = (n, m, 7), with the usual quantum 
numbers (n, /,m), while 7 labels the individual atoms in the unit cell. While the 
computation of the Fermi velocities is straight-forward e.g. in a plane-wave basis, their 
evaluation becomes tedious when using localized orbitals as required by many-body 
approaches such as LDA-I-DMFT. To this end, we employ the recently generalized Peierls 



approach [3l UHl [I9j , which we briefly summarize in Appendix A.l 

In their pioneering work, Rozenberg et al [20j analyzed the optical conductivity 
of V2O3 from the model perspective. It was concluded that the phenomenology of the 
temperature dependence in the conductivity can be understood by appealing to the 
physics of the one-band Hubbard model. In the current work, we will substantiate and 
extend these observations, based on a realistic multi-band setup. 

X In the case of a single orbital, vertex-corrections vanish in the limit of infinite lattice coordination [17j . 
Therewith, the response can be expressed in terms of spectral functions. In the general multi-orbital 
case this is an approximation. 
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3. Optical properties - results 



Our calculation of the optical conductivity is based on the previous LDA+DMFT 
electronic structure computation of Ref. |Tl], which used a one-particle Hamiltonian that 
was downfolded [211 122] to the vanadium t2g orbitals. Since optical experiments normally 
probe a much wider frequency range, we employ an upfolding scheme that includes 
higher energy states on the LDA level. Details are summarized in Appendix A. 2 




Figure 1. Theoretical optical conductivity of V2O3 at T « 390 K for a light 
polarization E \\ [a; z] = [0.13, 0.0, 0.041]. Contributions from different energy sectors 



(see [Appendix A. 2D : t2g t2g (red), t2g 



(green), 02p t2g (blue). 



Figure [T] shows our theoretical optical conductivity for V2O3 at T=390 K for 
the indicated light polarization. While in the Kohn-Sham spectrum, the t2g and 
bands are well separated, the correlations - accounted for by LDA+DMFT for the t2g 
orbitals only [H] - result in an intra-t29 conductivity that has weight up to energies well 
beyond the onset of transitions into the Cg at about 2 eV. Contributions stemming from 
transitions from the occupied oxygen 2p orbitals into the t2g occur from around 3.5 eV 
onwards. 

Before turning to a more detailed analysis of the different orbital contributions 
we compare our results to experimental data (see figure [2]). First, we notice the 
discrepancies between the experiments : Recent measurements on single crystals [23] 
agree well with previous single crystal experiments [20], but they are at variance 
with measurements using a polycrystalline film [23] . While the use of polycrystalline 
samples especially in a metal might be an issue, so is the fact that both single crystal 
experiments were performed up to energies of only a few eVs although the extraction 
of the conductivity involves a Kramers-Kronig transform. The low energy shape of 
the theoretical conductivity resembles the polycrystalline conductivity, but the absolute 
values differ. As to the single crystal one, we note that the order of magnitude compares 
favourably, while the shape tends to be comparable with the high temperature curves 
only. 
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At high energies, we see that both the onset and the shape of oxygen 2p derived 
contributions agree with experiment. The upfolding scheme that uses the 2p bands 
from LDA is thus a good approximation. This seems not to be true for transitions 
into orbitals : Compared with experiment |23], we reahze that spectral weight is too 
sharply defined and no identification of particular structures is possible. This calls for 
an LDA+DMFT calculation that includes all vanadium 3d orbitals. 

We now turn to a detailed analysis of orbital effects in the optical conductivity. 
This is a topic that has not at all been dealt with so far, since previous work neglected 
inter-band transitions altogether [23]. 

At low energy only a small Drude-like tail appears. This can be understood from 
the underlying electronic structure. Indeed, the metallic character of V2O3 is mainly 
a result from aig charge carriers that have spectral weight at the Fermi level only in a 
very limited region of the Brillouin zone (BZ), as can be seen in figure 4 of Ref. [llj. 

As can also be inferred from that work, the local spectral functions of aig and 
character display a pseudo-gap-like behaviour, and peak at finite energies rather than 
at the Fermi level, accounting for the feature seen at 0.5 eV in the conductivity. The 
latter orig ■uates from two types of transifoufl : At energies lower than 0.6 eV the 
spectral weight is mainly due to transitions from the aig into low lying orbitals, that 
are restricted to a small region in the BZ, whereas at slightly higher energies, 0.6 eV 
and above, contributions are in majority deriving from to transitions, which are 
possible in a wide region of the BZ, yet are less prominent at the F-point. 

At this point we again use our knowledge about the electronic structure of the 
compound: Reference [11] established an important orbital dependence of the quasi- 
particle coherence scale. Indeed, down to 390 K, excitations are far from being 
coherent : The imaginary part of the self-energy reaches —0.45 eV at the Fermi level, 
while aig excitations have reached their coherence regime in our calculation [TTl [3]. 
Therewith (ep aig carriers are (not) particularly sensitive to changes in temperature. 
As discussed above, the low-energy (< 0.6 eV) optical response is determined by ctig-e^ 
transitions, while above 0.6 eV e^^e^ transitions become dominant. Given these two 
facts, one can - even without explicit calculations - make some predictions about the 
behaviour of the optical response when the temperature is raised: Upon heating, the 
purely e^-derived contributions will not change as much as will those that involve the 
aig orbitals, so that the low-energy response will be more sensitive than the weight 
beyond 0.6 eV. In particular, a broadening (and thus reduction in height) is expected 
for the very low energy part. This gives a natural explanation for the dip behaviour 
that is observed in the experiments when the temperature is raised above ~ 450 K (see 
figure 1 in Ref. [23] or our figure [2]) [Jj]. Explicit calculations as a function of temperature 

§ The following is inferred from "momentum resolved optics" , i.e. from distinguishing contributions of 
different points in the BriouUin zone (not shown, see [3]). 

II Ref. [24j has invoked the change of the lattice constants upon heating in order to explain this dip. It 
is clear that the mechanism based on the orbital-selective coherence that emerges from our work could 
not have been observed in the calculation of Baldassarre et al, since, there, inter-band transitions were 
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(including inter-band transitions) would be desirable to confirm the picture emerging 
from our results. This challenging project is left for future work. 




1 2 3 4 5 6 

co[eV] 

Figure 2. Theoretical and experimental conductivity of V2O3. Theory (same 
polarization as in figure [T]) : total (bold black), orbital contributions (shades of gray, 
dotted). Experiments : polycrystalline film, (T=200, 400 K) Qazilbash et al [23] (blue, 
green, long dashed), single crystal (T=300-600 K) Baldasarre et al [M] (blue to red, 
short dashed). 
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Figure 3. Theoretical optical conductivity of (Vi_a;Cr2;)203, x=3.8% at T = 580 K 
for a light polarization E \\ [a; z] = [0.13, 0.0, 0.0415]. Contributions from different 
energy sectors : t2g — > t2g (red), t2g (green), 02p <2g (blue). Gray lines show 

the results of pure V2O3 (figure [T]) for comparison. 

In figure [3] we show theoretical results for the insulator (Vi_^Cra;)203 (x=3.8%). As 
discussed previously [iTj, we have low but finite spectral weight at the Fermi level, which 
results in some optical weight at low frequency. Unfortunately, no experimental data 
is available for this composition. Compared with pure V2O3, we note the suppression 

neglected (see their footnote 26 and our discussion in [Appendix A.l[ ). 
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of low energy spectral weight and the clear distinction of transitions into the t2g upper 
Hubbard bands at ~4 eV. 



4. Conclusions 



In conclusion, we presented calculations of the optical conductivity of V2O3 in 
its paramagnetic phases using the generalized Peierls approach. We obtain good 
agreement with experiment and propose an explanation for the experimentally witnessed 
temperature dependence of the response as a signature of the orbital selective coherence 
of the system. Our upfolding scheme to include higher energy orbitals captures well the 
transitions involving oxygen states, but reveal the necessity to include all 3d orbitals in 
an LDA+DMFT electronic structure calculation for this compound. 
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Appendix A. Details of the formalism 

Appendix A.l. Generalized Peierls substitution approach to Fermi velocities 
The Fermi velocities in ([1]) are given by elements of the momentum operator V : 

v^t = -{\iL'\V^\kL) (A.l) 
m 

L={n,l,m,'y) and 7 labels atoms in the unit-cell. In plane waves flA.l|) is easily 



calculated, while using a localized Wannier-like basis XRL(r) = (r|RL)=^,^ e~*''^(r|kL) 
renders the evaluation tedious. Inspired by the Peierls substitution approach [16] of 
lattice models, we can separate the above into [HI [THl [19] : 

C = I {9,H^' - <Pl' - Pl)H^'") + -^H [{XrJ] (A.2) 

The terms in brackets are the Fermi velocity in the Peierls approximation, which is here 
generalized to multi-atomic unit-cell : pi is the a-component of the position of atom 
7 within the unit-cell. This velocity is easy to evaluate since it involves only elements 
of the Hamiltonian. Crucial is that ^jj[{Xrl}] reduces to intra-atomic contributions 
in the limit of strongly localized orbitals Xni [2]' which makes the generalized Peierls 
velocity a good approximation for e.g. 3d and 4f systems. 

Finally, we ask if the computation of the Fermi velocities is really necessary 
in practice, or if one could also resort to a simpler approximation consisting in 



Multi- orbital Effects in Optical Properties of Vanadium Sesquioxide 



7 



simply omitting the Fermi velocities. Due to its simplicity, this approximation is in 
fact relatively popular to obtain qualitative trends of optical properties in correlated 
systems [25l [2^ . Since then the conductivity is then a simple convolution of spectral 
functions, inter-band transitions {L^L') are neglected and intra-band transitions not 
properly weighted. As an illustration, we show in figure lAll a comparison of the optical 
conductivity calculated within the generalized Peierls formalism compared to the one 
computed from the simple convolution of spectral function^: Besides the obvious 
discrepancy in absolute value, omitting the Fermi velocities results in a noticeable change 
in shape, too. This is owing to the momentum dependence of the matrix elements that 
favours certain regions in the Brioullin zone, while attenuating others. 



T=390K, E||[xOzl 



generalized Peierls, total 
Vt,„^Vt„ 




1 2 3 4 5 6 7 
a)[eV] 



Figure Al. Generalized Peierls conductivity versus "no velocity approximation". 



Appendix A. 2. Up folding scheme for higher energy transitions 

Although the treatment of many-body correlations can often be cast into an effective 
low energy, downfolded system, the range of its validity is usually far exceeded by 
optical measurements. Thus it is desirable to allow for optical transitions into higher 
energy orbitals. Also, the computation of the Fermi velocities and the downfolding 
procedure do not commute [TB], and hence it makes a difference to which Hamiltonian 
the generalized Peierls approach is applied. As a matter of fact, the Wannier functions of 
a full Hamiltonian are more localized than for the downfolded one, whereby the Peierls 
approximation becomes more accurat^. The key quantity for the conductivity is the 
orbital trace in ([1]). For any unitary transformation ?7k holds 

ti{v^A^{u')v^A^{u'+u)} = ti {UlvkUi,Ak{u')Ulvi,UkAi,{u' + u)} (A.3) 

where we defined A]^ = f/^^kf^;- In the case of a pure band-structure calculation (no 
self-energy), we can choose the transformation such that it performs the downfolding, i.e. 
the spectral functions Aif_ acquire a block-diagonal form. We shall distinguish between 

^ In order to have comparable scales, we chose for the latter case Wk = tqI, with the Bohr-radius vq. 
+ Indeed the downfolding can be viewed as a unitary transformation that block-diagonalizes the 
Hamiltonian (see below). The change in accuracy manifests itself in the basis dependence of the optical 
conductivity within the Peierls approach. 
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the low energy (=Sf) and the high energy block (t^) : An LDA+DMFT calculation 
will add local Coulomb interactions only to the former after the block-diagonalization, 
which results in a self-energy that lives in this sub-block, while high energy bands remain 
unchanged and the block-diagonality is retained. Clearly the downfolding procedure is 
not exact in the many-body framework. Indeed the matrices that block-diagonalize the 
true interacting system also depend on frequency due to the dynamical nature of the 
self-energy. Yet, when granting the approximative validity of the downfolding, and using 
the [/k of the band-structure calculation, we can specify 






with = t^k'^kf^k- The many-body spectra L, L are substituted into the =Sf -sector, 
while H, H of the ^^-sector stem from the initial band- structure, and (]A.3|) reads 

LViLVi + LWHW^ + HV2HV2 + HW^LW (A.5) 

For transitions within the =Sf-block, the velocity Vi appears, which is the ^-block of the 
transformed velocity. It is different from the element computed after the downfolding. 
With the above, we moreover include transitions from, to and within the high energy 
bloclfl. Comparing to experiments then allows to assess whether correlation effects 
substantially modify also the spectrum of downfolded orbitals, or whether for them the 
initial band-structure is satisfying (see above for the V2O3 case). 
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